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Abstract. The gravitational collapse of a star is an important issue both for general relativity and 
astrophysics, which is related to the well known "frozen star" paradox. Following the seminal work 
of Oppenheimer and Schneider (1939), we present the exact solution for two dust shells collapsing 
towards a pre-existing black hole. We find that the inner region of the shell is influenced by the 
property of the shell, which is contrary to the result in Newtonian theory and and the clock inside 
the shell becomes slower as the shell collapses towards the pre-existing black hole. This result in 
principle may be tested experimentally if a beam of light travels across the shell. We conclude that 
the concept of the "frozen star" should be abandoned, since matter can indeed cross a black hole's 
horizon according to the clock of an external observer. Since matter will not accumulate around the 
event horizon of a black hole, we predict that only gravitational wave radiation can be produced in 
the final stage of the merging process of two coalescing black holes. Our results also indicate that 
for the clock of an external observer, matter, after crossing the event horizon, will never arrive at the 
"singularity" (i.e. the exact center of the black hole). 
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INTRODUCTION 

The "frozen star" paradox is a well known novel phenomenon predicted by general 
relativity, i.e. a distant observer (O) sees a test particle falling towards a black hole 
moving slower and slower, becoming darker and darker, and is eventually frozen near 
the event horizon of the black hole. However, as we discussed in [1], several possible 
explanations to this phenomenon have been proposed, but none of these is completely 
satisfactory. In the previous work, we found that an external observer should be able to 
observe matter falling in a black hole, based on our stationary solution of a dust shell 
around a pre-existing black hole. In this paper we further obtain the exact dynamical 
solution for two shells collapsing towards a pre-existing black hole, in particular to find 
the difference between this solution and that for the gravitational collapse of a uniform 
dust ball and the implication for the "frozen star" paradox. 
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THE DYNAMICAL SOLUTION 



We follow the method adopted in the seminal work of Oppenheimer and Schneider 
(1939) [211], i.e. obtaining the solution in the comoving coordinates first and then trans- 
forming it into the ordinary coordinates. 

The general form of the metric in the comoving coordinates is (G = c = 1 is adopted 
throughout this paper) 

ds^ = dx^ - e^dR^ ~e''{de^ + sin^ ed(l)^), (1) 

where G){R, t) and Co{R, t) are the unknown functions. In the comoving coordinates, the 
only non-zero component of the energy momentum tensor is = p. As the results in 
[0], we can obtain the following equations from the field equation 

e'^ = (FT + G)4/3, (2) 

= e"w'V4, (3) 

Snp = ^ {t + G/F)-'{t + G'/F')-\ (4) 

where F and G are arbitrary functions of R. Dot and prime are the partial differentiation 
with respect to t and R, respectively. At t = 0, we obtain 

F = -^6;rpo(G2 + Ci), (5) 

where Ci is the integration constant. As in ['2'], we choose G = R^/^. 

As obtained in |l2(], the value of F in region V is —jy/ro, where ro = 2M (M is the 
total mass of the system). Then to assure the solutions are continuous at the boundaries, 
the results of F in other regions should be 




-l^J^{R^-al)+ro ;4<i?<a2 

where ro = 2M, = 2m, /q =2{m-\-m\). 

Up to now we have obtained the solution in the comoving coordinates. Then we try to 
transform the solution into the ordinary coordinates, in which the metric has the form as 

ds^ = B{r, t)dt^ -A{r,t)dr^ - r^{de'^ + sin^ ed<^'^). (7) 

The transformation of r is obvious. We must choose r = e^l^ = {Ft + G)^/^. 
Using the contravarient form of the metric and requiring that the g^' term vanishes, 
we have 

t'/i = rr. (8) 




FIGURE 1. (a) is the configuration of the dynamical solution for two shells in comoving coordinates. 
m, mi and m2 are the masses of the black hole, the shell 1 and 2, respectively. a\ and oi are the radii of the 
inner and outer boundaries of shell 1, respectively, and 02 are the radii of the inner and outer boundaries 
of shell 2, respectively, (b) is evolution curves of the dynamical solutions for 02 — lOro, 02 = 8ro, oi = 5ro, 
a\ — 2.5ro, '"0 = 1 /3ro, and r'^ = 2/3ro with coordinate time, (c) is the comparison of the evolution of 
the outer boundary of shell 1 between the case with (solid) or without (dashed) shell 2. The parameters 
of the shells are 02 = lOro, = 6ro, a\ = 2.5ro, fli = 5ro, = l/5ro, and = 2/5ro, where ro is the 
Schwarzschild radius corresponding to the total mass of the system, i.e., the sum of the masses of the 
pre-existing black hole, shell 1 and shell 2. The inset is the ratio of the solid line to the dashed line. It is 
seen clearly that outside matter, though spherically symmetric, does influence the motion of matter inside. 



Substituting the expressions of r in five regions into Eq. ([8]), we can obtain five partial 
differential equations of t{R, t). 

In region V, the solution of Eq. ([8]) is 

t = LW, X = -^(i?3/2 _ ,3/2) _ 20^ + roln ^^±^, (9) 

where L is an arbitrary function of x. Since the metric in the ordinary coordinates in 
region V is the Schwarzschild metric, L{x) should be x. 

In general case, the partial differential equations in other regions can be solved 
numerically only, using the solution in region V as the boundary condition. 

After obtain the transformations in every region, we could know how the shells evolve 
with the coordinate time. Fig. la and b show the configuration of this problem and the 
evolution of the shells in the ordinary coordinates. In Fig. Ic, we show the evolution 
curves of the outer boundary of shell 1 in the case with or without shell 2. The outer 
boundary of shell 1 will take more time to approach the asymptotic line in the case with 
shell 2. 

DISCUSSIONS AND CONCLUSIONS 

In the Newtonian gravitation theory, if the matter is spherically symmetric, the outer 
matter will not influence the inner region. However, as shown in Sec. II, in general 
relativity case, the clock in the inner region of the shell is slower compared with the case 
without the shell. This effect may be testable experimentally in principle, e.g. if a beam 
of light travels across the shell and a parallel beam of light passes outside the shell (far 



enough from the shell), then the two beams of light will travel at different velocities with 
respect to the observers outside the shell. 

The inner shell can cross the Schwarzschild radius in the two shell case. In this sense 
we could observe the matter falls into a black hole and the "frozen star" paradox dis- 
cussed in [Sl] is naturally solved. As pointed out previously by us, the origin of the 
"frozen star" is the "test particle" assumption for the infalling matter, which neglects 
the influence of the infalling matter itself on the metric [1]. In fact, for some practical 
astrophysical settings, the time taken for matter falling into a black hole is quite short, 
even for the external observer We find that the infalling matter will not accumulate 
outside the event horizon, and thus the quantum radiation and Gamma Ray bursts pre- 
dicted in [4] and [5] are not likely to be generated. We predict that only gravitational 
wave radiation can be produced in the final stage of the merging process of two coalesc- 
ing black holes. Future simultaneous observations by X-ray telescopes and gravitational 
wave telescopes shall be able to verify our prediction. 

It is also interesting to note that, as can be seen from Fig.IIh, in ordinary coordinates, 
the matter will not collapse to the singularity {r = 0) even with infinite coordinate time 
(if Tq = 0, the inner boundary of the inner shell will take infinite time to arrive at r = 0). 
It means that in real astrophysics sense, matter can never arrive at the singularity (i.e., 
the exact center of the black hole) with respect to the clock of the external observer. 
Therefore, no gravitational singularity exists physically, even within the framework of 
the classical general relativity. 
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